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ABSTRACT 

Context. Four clusters of red supergiants have been discovered in a region of the Milky Way close to base of the Scutum-Crux Arm 
and the tip of the Long Bar Population synthesis models indicate that they must be very massive to harbour so many supergiants. If 
the clusters are physically connected, this Scutum Complex would be the largest and most massive star-forming region ever identified 
in the Milky Way 

Aims. The spatial extent of one of these clusters, RSGC3, has not been investigated. In this paper we explore the possibility that a 
population of red supergiants could be located in its vicinity. 

Methods. We utilised 2MASS JHK S photometry to identify candidate obscured luminous red stars in the vicinity of RSGC3. We 
observed a sample of candidates with the TWIN spectrograph on the 3.5-m telescope at Calar Alto, obtaining intermediate-resolution 
spectroscopy in the 8000-9000A range. We re-evaluated a number of classification criteria proposed in the literature for this spectral 
range and found that we could use our spectra to derive spectral types and luminosity classes. 

Results. We measured the radial velocity of five members of RSGC3, finding velocities similar to the average for members of 
Stephenson 2. Among the candidates observed outside the cluster, our spectra revealed eight M-type supergiants at distances < 18' 
from the centre of RSGC3, distributed in two clumps. The southern clump is most likely another cluster of red supergiants, with 
reddening and age identical to RSGC3. From 2MASS photometry, we identified four likely supergiant members of the cluster in 
addition to the five spectroscopically observed. The northern clump may be a small cluster with similar parameters. Photometric 
analysis of the area around RSGC3 suggests the presence of a large (> 30) population of red supergiants with similar colours. 
Conclusions. Our data suggest that the massive cluster RSGC3 is surrounded by an extended association, which may be very massive 
(k, 10 5 M©). We also show that supergiants in the Scutum Complex may be characterised via a combination of 2MASS photometry 
and intermediate-to-high-resolution spectroscopy in the Z band. 

Key words, stars: evolution - early type - supergiant - Galaxy: structure - open clusters and associations: individual: Alicante 7 - 
RSGC3 - Stephenson 2 
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^ ^ r , , , the most heavily obscured, and also the young est at an estimate d 

quence RSGs are powerful beacons that can be identified even T=u±2 M and M 3±1 x 10 4 M (Dav i ^t^i] [2008). 

behind huge amounts of obscuration. Unfortunat ely, the physical Alicante g located in it s immed iate vicinity, se ems slightly 

properties of RSGs are still poorly constrained jLeyg sque et alj older and father smaUer dNegueruelaet al.ll201toll . RSGC2 = 

|200l), although they are known to span a wide range of lu- St henson 2 (Ste 2) is the least obscured and apparent i y most 

minosities, depending on their initial masse s and evolutionary massive of a ll with T=17±3 Myr and Minitial =4±l x 1O 4 M 

S' t'u P ym / lo | (Lbol/Le) ~ 4 -°- 5 - 8 <Meynet&Maedej toavi e S et alJ l200i henceforth fe. Finally, RSGC3 lies at 

pi Their infrared spectra, however, do not reflect this va- somg distance from aU ^ Qthers and has a n estimated r=16- 

nety. A 1 supergiants with spectral types between mid G and late 20 Myr and an in ferred M initial =2-4xl0 4 M o dClark et al.l2009at 

M display very similar infrared spectra, with a strong degen- | A lexander et al.l l2009). Collectively, these clusters host a signif- 

eracy between temperature and luminosity. As a consequence, ic£fflt fraction of ^ knQwn RSG lation in the Gal 

absolute magnitudes cannot be interred from mirared spectra, _ . . , 1T ^o^i ^o^^x 

implying that it is not possible to assign physically meaningful Th , e ^ of sl S ht t0 the ^ SG c J usters ^ RSGCs > Passes through 

parameters to any individual obscured RSG. ^ several Galactic arms and reaches an area of very high ob- 

In the past few years, a flurry of discoveries has revealed sev- scuration at a distance comparable to those of the clusters 

eral clusters of RSGs located in a small region of the Galactic (Negueruela^|2010a|). Stellar densities are extremely high at 

moderate K magnitudes, meaning that the cluster populations are 

* Partially based on observations collected at the 3.5-m telescope very hard to identify, except for the RSGs. Because of this, mem- 

(Calar Alto, Spain) bership of individual stars is difficult to assess and the actual ex- 
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tent of the clusters is uncertain. Membership in RSGC1 and Ste 2 
has been defined in terms of a common radial velocity, v ra d, mea- 
sured from high-resolution K-ba.nd spectra (iDavies et al .1 12007, 
2008). In the case of RSGC1, the obscuration is very high and 
the RSG members are intrinsically very bright, resulting in a rel- 
atively easy selection of candidate members. On the other hand, 
the line of sight to Ste 2 is very rich in bright red stars. lD07l 
find a large number of foreground giants and also several pu- 
tative RSGs with v ra d apparently excluding membership, even 
though their analysis was confined to a circle with radius r = 7' 
around the nominal cluster centre. For the other clusters, dynam- 
ical analysis has not yet been carried out and the determination 
of membership is less secure. 

With the data currently available, the exact nature of this 
Scutum Complex is unclear. There have been suggestions that it 
represents a giant star formation region triggered by the dynam- 
ical excitation of the Galactic Bar, whose tip is believed to inter- 
se ct the Scutum-Crux Arm close to this region (see discussion 
in lD07l : iGarzon et al.lfl997l) . If all the RSG clusters are located 
at a common nominal distance of 6.5 kpc, they span ~ 500 pc in 
projection. Alternativel y, if there is a gian t Molecular Ring at the 
end of the Bar (e.g., Rathbor ne et al .1 12009), our sightline could 
cut its cross section over a distance ~ 3 kpc. 

Heretofore, all studies of the RSGCs have been carried out 
solely with infrared data. Intermediate-resolution K-ba.nd spec- 
troscopy has been used to observe the CO bandhead absorp- 
tion features near 2.3/zm, which identify cool, evolved stars. The 
equivalent width (EW) of the first bandhead feature is a valid 
estimator of spectral type at intermediate resolution, though the 
relationship is bimodal (showing a similar slope for giants and 
supergiants, but higher EWs for the latter), a n d scatter is large 
(Da vies et al] 120071; iNegueruela et al.l l2010al) . I Alexander et al] 
( 2009) explore a similar relation for low-resolution spectra. In 
both cases, the accuracy in the spectral-type determination is 
estimated at ±2 subtypes. Until the advent of infrared spectro- 
graphs with multiplexing capabilities, extending this methodol- 
ogy to large samples will consume large amounts of telescope 
time. 

Optical spectrographs with high multiplexing factors do, 
however, exist. Highly obscured RSGs are, unfortunately, too 
faint to be detected in the optical range. The far red optical re- 
gion, between 8000 A and 10000 A offers a good compromise, 
as it can be reached with optical instrumentation and has lower 
extinction. B y using existing ca talogues, namely B/DENIS and 
USNO-B 1 .0 (IMonet et al.l2003l) . we have verified that more than 
50% of the confirmed RSG members of Ste 2 and RSGC3 are 
brighter than i - 16 and can be observed at intermediate-high 
resolution with existing instrumentation. 

In this paper, we present a pilot study to address several is- 
sues that will allow us to evaluate the possibility of detecting 
RSGs over a wide area by using optical spectrographs. Firstly, 
we explore the use of photometric criteria and catalogue cross 
correlation to detect candidate luminous red stars. Secondly, we 
re-evaluate several criteria presented in the literature to estimate 
spectral types and luminosities with intermediate-resolution far 
red optical spectra. Finally, we assess the existence of an 
extended population of RSGs around RSGC3 that may jus- 
tify wide-field searches. The paper is organised as follows. In 
Section O we set photometric criteria for selection of candi- 
dates and describe the observations of a sample of candidates. 
In Section^ we discuss several methods of classification for red 
stars in the near-IR. In Section HI we present the analysis of our 
observations. Finally, we discuss their implications in Section [5] 



2. Target selection and observations 

RSGC3 was discovered as a concentration of very b ri ght infrared 
sources (IClark et a l. 2009a; Alexander et al. 2009). Clark eiaL| 
( 2009a|) explored a circle of radius r = 7' around the nominal 
centre of RSGC3, finding a strong concentration of RSGs to- 
wards the centre, but also several photo metric candidate mem- 
bers at the edges of their search area. lAlexander et al.l (2009) 
analysed a smaller area. 

In regions of high extinction, traditional photometric criteria 
cannot always be employed to select stars of a particular spec- 
tral type and luminosity class. While RSGs are very bright in 
K, several classes of red giants are also bright infrared sources, 
but much more numerous. Despite this, if we are to detect RSGs 
in the area surrounding RSGC3, we might expect them to have 
similar photometric characteristics to known cluster members. 

2.1. Photometric criteria 

We have used 2MASS JHK S photometry dSkrutskie et al.ll2006l) 
to ident ify possible RSG stars. As discussed in Negueruela et al. 
(I2010al) . the reddening-free parameter O = (J-H)-l.8x(H-K s ) 
allows the separation of early-type stars (Q ^ 0.0) and late-type 
main-sequence stars and red clump giants (Q ^ 0.4). RSGs show 
a peculiar behaviour and present very different Q values. Most of 
them deredden to the position of yellow stars, with Q « 0.2-0.3, 
though a fraction fall close to the position of red stars (Q « 0.4). 

For this pilot programme, we selected stars within 20' of the 
nominal centre of RSGC3, with (/ - ^s) > 1.3. The colour cut 
was chosen to be redder than the intrinsic colours of any nor- 
mal, unreddened star, as we expected our objects to suffer sig- 
nificant extinction. On the other hand, we did not want to force 
candidates to have the same colours as cluster members, since 
extinction might be variable on small scales over this field. We 
selected stars brighter than = 7.0, as RSGs are expected to 
be bright. Typical magnitudes for confirmed members of RSGC3 
lie in the = 5 - 6 range, and so = 7.0 allows for much 
heavier extinction. Finally, we imposed 0.1 < Q < 0.4, the range 
where most RSGs lie, in order to exclude foreground red-clump 
giants and main- sequence M-type stars. We note that the criteria 
are unlikely to select all RSGs, as some of them have colours 
Q « 0.4, or may be so heavily reddened as to have > 7.0. 
The criteria were chosen to select most RSGs, while leaving out 
many likely interlopers. 

We then checked our selection against the B/DENIS and 
USNO-B 1.0 catalogues. Stars clearly detected in the DSS2 blue 
plates were eliminated from the sample. In addition, we set a cut 
at i > 9.5, as we did not expect our targets to be brighter than 
this limit, given that extinction is much higher in i than in K. The 
size of the telescope forces us to a limit i < 15.5. This limit was 
never reached during the observations, as the seeing was rela- 
tively poor throughout the run. Only stars brighter than i = 14.5 
were observed. 

There are ~ 50 stars fulfilling these criteria and lying at dis- 
tances > 7' and < 20' from the cluster centre. They are not 
evenly spread, and we can identify three areas of high density, 
which we select for systematic observations (see Fig. [TJ. Field 1 
lies 14' South- West of RSGC3. There is a compact group of can- 
didates with several other objects around. Field 2 contains an 
obvious concentration of very bright (in K) stars, most of which 
pass all the photometric criteria. It is located ~ 12' North- West 
of the centre of RSGC3. We also include a few fainter stars a 
few arcmin to the North. Finally, Field 3 does not appear to con- 
tain any obvious aggregate, instead comprising an area of uni- 
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Fig. 1. Colour-composite 2MASS JHK$ finder of the area surveyed. The image covers ~ 55' x 35' (roughly corresponding to 
100x60 pc at 6 kpc). The circle showing the position of RSGC3 has radius 6'. Yellow circles show stars that pass all the photometric 
criteria in S ection l2.1l Red circles indicate stars which, in addition, are detected as DENIS or USNO.B-1 sources with i> 13.2 (see 
Section l53l for this selection criterion). Crosses indicate the stars spectroscopically identified as reddened RSGs in this paper. North 
is up and East is left. 



formly distributed photometric candidates. It lies -12' South- 
East of the cluster core. All the photometric candidates are listed 
in Table [U where they are given an identifier of the form FnSnn, 
where the first digit identifies the field and the last two digits 
identify the individual star. 



2.2. Spectroscopy 

Observations were carried out with the Cassegrain Twin 
Spectrograph (TWIN) of the 3.5 m telescope at Calar Alto 
(Almeria, Spain), during a run on 2009 July 10-11. We used 
the red arm equipped with the SITe#20b_12 CCD and the T06 
(1200 In mm" 1 ) grating. The grating angle was chosen to cover 
the 8000-9000 A range. The 1'/ 2 slit gives a resolution element 
of 1.2A, measured on the arc lamps. This is equivalent to a re- 
solving power R = 7 000. For this pilot run, we chose a longslit 
spectrograph, as it allows a much better selection of individual 
targets and observation of stars in the crowded cluster cores. 

In addition to the stars in the three fields listed abov e, we 
observed seven objects from the list of members of Ste 2 in D07 
to be used as references for radial velocity measurements, and 
five photometric members of RSGC3. 

The seeing was variable, with typical values ~ 175 (the 
monitor was not active at the time, so only rough estimates 
were available). Transparency was high, except for the first two 
hours of the first night, when some clouds were present close 



to the horizon. Typical exposure times ranged from 150 s for 
the brightest objects (i ^ 10) to 2000 s for objects fainter than 
i = 14.0. To reduce integration times, most targets were observed 
in pairs, aligned along the slit. Because of this, parallactic an- 
gles were not used. ThAr arc lamps were observed immediately 
before or after any exposure, and each time that the telescope 
was moved more than a few arcmin. All exposures have been 
calibrated with, at least, one arc taken less than 1800 s before 
or after the midpoint of the observation and, in many cases, are 
sandwiched between two arcs. 

All the data have bee n reduced using the Starlink 
software packages ccd pack (iDraper et al.l l2000b and figaro 
(Shortridge et al.l 1 19971) and analysed using figaro and dipso 
( Howarth et al ■1 11998b . Arc frames were extracted at the position 
of each target star, and the respective solutions checked. Typical 
RMS values for the fit to a third-degree polynomial are « 0.1 A. 
Solutions at different positions within a given frame are signifi- 
cantly shifted, with differences up to two pixels. 

The average signal-to-noise ratio (SNR) of our targets was 
not high, typically 50-70 per pixel, though higher for the bright- 
est targets, and lower for some objects (a few spectra have 
SNR- 35). As an example, in Fig. |2l we show the spectrum 
of star D18 (i.e., star #18 from|D07|) in Ste 2, which has a SNR 
~ 90 per pixel (~ 145 per resolution element), with a number of 
featu res identified, after different references i n the literature, no- 
tably iKirkpatrick et al.l dl99lb . lMaliikl (Il997h . and lCenarro et al.l 
( 2001J - " 
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Table 1. Candidates selected as possible red luminous stars, with their characteristics. 



ID 


2MASS 


J 1 




^s 1 


Q 




Spectral 


Spectral 


E(J-Ks) 


Vlsr 


















Type (TiO) 


Type 




(kms" 1 ) 


Field 1 


F1S01 


18445403- 


-0329009 


7.55+0.02 


6.26+0.04 


5.71+0.02 


0.30 


11.08+0.02 


M7 


M7 Ib/II 


0.55 


+43 


F1S02 


18444686- 


-0331074 


8.78±0.03 


6.71+0.02 


5.76+0.02 


0.38 


14.44+0.04 


Ml 


Ml la* 


1.96 


+73 


F1S03 


18444439- 


-0334197 


8.21+0.02 


6.86+0.04 


6.32+0.02 


0.39 


11.30+0.02 










F1S04 


18443941- 


-0330003 


8.66+0.03 


6.76+0.04 


5.82+0.02 


0.19 


13.99+ 0.03 


M2 


M2Iab* 


1.74 


+ 112 


F1S05 


18443100- 


-0330499 


9.87+0.02 


7.93+0.05 


6.97+0.03 


0.21 


15.44+0.06 










F1S06 


18442945- 


-0330024 


8.45+0.03 


6.39+0.03 


5.41+0.02 


0.30 


14.13+0.04 


M3 


M2Ia* 


1.94 


+93 


F1S07 


18442796- 


-0329425 


8.06+0.02 


5.91+0.04 


4.85+0.02 


0.25 


14.46+0.04 


M5 


M4Iab* 


2.04 


+93 


F1S08 


18443037- 


-0328470 


7.34+0.02 


5.20+0.04 


4.10+0.30 


0.16 


13.70+0.03 


M4 


M4Ib* 


2.07 


+75 


F1S09 


18442053- 


-0328446 


9.28+0.02 


7.24+0.07 


6.21+0.03 


0.21 


15.32+0.06 










F1S10 


18442322- 


-0324578 


8.08+0.03 


6.64+0.05 


5.94+0.02 


0.18 


11.91+0.02 


M6 


M5 Ib/II 


0.95 


+59 


fish 


18441700- 


-0325360 


8.68+0.03 


7.40+0.04 


6.89+0.02 


0.36 


11.44+0.02 










F1S12 


18444065- 


-0335199 


8.90+0.03 


7.49+0.03 


6.90+0.02 


0.37 


12.22+0.03 










Field 2 


F2S01 


18444846- 


-0313495 


9.14+0.02 


7.33+0.03 


6.50+0.02 


0.32 


14.01+0.04 


Ml 


M0 la* 


1.63 


+72 


F2S02 


18444431- 


-0313350 


8.69+0.02 


6.39+0.03 


5.23+0.02 


0.22 


15.52+0.07 










F2S03 


18444044- 


-0314142 


8.06+0.02 


6.22+0.03 


5.31+0.02 


0.21 


13.21+0.03 


Ml 


Ml lab* 


1.67 


+86 


F2S04 


18444283- 


-0315168 


8.53+0.02 


6.69+0.04 


5.80+0.02 


0.24 


13.69+0.03 


M2 


Ml. 5 lab* 


1.64 


+87 


F2S05 


18444023- 


-0315329 


7.21+0.03 


5.51+0.03 


4.72+0.03 


0.30 


12.40+0.03 


M8 


M7 II 


1.19 


+61 


F2S06 


18442823- 


-0312558 


7.50+0.03 


6.38+0.03 


5.94+0.02 


0.34 


9.70+0.03 










F2S07 


18443833- 


-0316596 


7.32+0.02 


6.05+0.04 


5.49+0.02 


0.26 


10.68+0.03 


M7 


M6III 


0.58 


+85 


F2S08 


18443324- 


-0306334 


8.44+0.03 


7.34+0.06 


6.88+0.04 


0.27 


10.73+0.02 










F2S09 


18443220- 


-0304597 


9.04+0.03 


7.29+0.05 


6.44+0.04 


0.23 


13.77+0.03 


M5 


M5.5III 


1.36 


+94 


F2S10 


18443172- 


-0302217 


8.40+0.02 


6.85+0.03 


6.19+0.03 


0.35 


12.15+0.03 










F2S11 


18444658- 


-0302416 


8.65+0.02 


6.69+0.03 


5.79+0.02 


0.33 


14.18+0.04 










Field 3 


F3S01 


18461105- 


-0330091 


7.08+0.02 


5.79+0.03 


5.26 +0.02 


0.33 


10.04+0.04 


M3 


M4II-III 


0.65 


-22 


F3S02 


18460994- 


-0334581 


8.10+0.02 


6.80+0.05 


6.26+0.02 


0.33 


10.94+0.02 


M2 


M3.5III 


0.70 


+ 10 


F3S03 


18455541- 


-0331254 


8.37+0.02 


6.82+0.05 


6.17 +0.02 


0.38 


13.04+0.02 


M8 


M7III 


0.90 


+66 


F3S04 


18455357- 


-0332069 


9.30+0.02 


7.31+0.04 


6.38+0.03 


0.31 


14.76+0.04 










F3S05 


18454257- 


-0332406 


6.68+0.02 


5.27+0.02 


4.60+0.02 


0.19 


11.67+0.02 


M9 


M7III 


0.78 


-37 


F3S06 


18453198- 


-0331536 


7.99+0.02 


6.31+0.03 


5.51+0.03 


0.24 


11.5+0.3 










F3S07 


18453212- 


-0334044 


8.40+0.03 


6.45+0.03 


5.51+0.03 


0.25 


13.4+0.3 










F3S08 


18453872- 


-0336069 


7.55+0.03 


6.23+0.033 


5.66+0.02 


0.29 


10.96+0.02 


M6 


M5 II-III 


0.68 


+59 


F3S09 


18453660- 


-0339191 


7.20+0.02 


5.83+0.04 


5.24+0.02 


0.30 


10.42+0.03 


M3 


M4II 


0.79 


+35 


F3S10 


18460233- 


-0338568 


8.54+0.03 


6.10+0.04 


4.94+0.02 


0.36 


15.31+0.05 










(1) JHK S 


magnitudes from 2MASS. i magnitudes from DENIS, except for F3S06 and F3S07, which are photographic / from USNO. 





Likely members of the RSGC3 association. 



3. Spectral classification 

The classification of red luminous stars in the far red region has 
been explored by several authors. Different criteria have been 
proposed dependi ng on the resolution av ailable. For example, 
Bos chi etalT(l2QQ3l), using echelle spectra (iMunari & Tomasellal 
1999; Marrese et al.H 2003). propose classification criteria based 
on ind ivid ual line ratios. At low er resolution, Gin estet et al.l 
(11994 and lCarquillat et al.l (1 19971) investigated the variation of 
different line strengths with temperature and luminosity with the 
intention of finding valid criteria for the determination of effec- 
tive temperature (T e ff) and luminosity. 

In general, at a given metallicity, the st rength of the metal- 
lic spectrum increas es with luminosity (iGinestet et al.l Q994; 
ICarquillat et al.|[l997l) . An important case in point is the strength 
of the Can triplet. These lines dominate the /-band spectra of 
stars with spectral types between early F a nd M3 and their 
strength is s trongly correlated with luminosity (iDiaz et al.f i 9 89: 
ZhouT ll99ll) . although they also demonstrate a dependence on 
metallicity and spectral type. Several authors have studied the 
dependence of the equivalent width (EW) of these lines with 



luminosity class, but the results depend to some degree on the 
resolution used a nd the way in which the continuum is defined 
rcf. lMalliklfT997h . 

While main sequence stars are easy to tell apart, differences 
between giants and supergiants are more subtl e. Amongst sev- 
eral line strengths and line ratios studied by Carqu iilat et al.l 
(11997b . the strength of the Tii and Fei blend at 8468A seems 
to be the best discriminant at a given spectral type. However, all 
luminosity criteria show a certain degree of dispersion due to 
differences in chemical composition. In addition, many metal- 
lic lines show dependence on both temperature and luminosity. 
As a consequence, the bo undary between luminosity classes I I 
and lb may be uncertain (Gineste fet al.|[l9 94; Zhu et al.lll999l) . 
in particular for late K types (IGinestet et al.ll 1994b . 

The different luminosity calibrations in the literature have 
been derived from spectra with different resolutions, using di- 
verse measuring methods. Moreover, all these works have used 
spectra with high SNR. As our spectra have moderately high 
resolution and only moderate SNR, we choose to measure the 
EW of lines relative to the local continuum, using windows de- 
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00 00 

I 



00 



II II 




Wavelength (A) 

Fig. 2. Intermediate resolution spectrum of star D18 in Ste 2, 
with features identified, and the position of the TiO bandheads 
used for spectral classification of later- type stars marked. This 
object is Ml. 5 lab according to the features in this spectral re- 
gion. Also shown (displaced upwards by 0.4 cont inuum units) 
is the sky transmission in this spectral region, from Hinkle et al. 
(20031). 



fined by the apparent line width (which, for example, turn out to 
be ^14A for Can 8498 A and ^18A for Can 8542A in all our 
stars). As there is an amount of subjectivity in this definition, 
we estimate the internal errors in the measurements by taking 
a large number of measurements with small differences in the 
definition of the continuum and integration limits. Typically the 
errors are found to be « 3%. T o calibrate our mea surements, 
we utilised the spectral library of Cenarr o et al .1 (120011) since the 
spectra in this library have similar resolution to oursQ. In addi- 
tion, many of the stars included are defined as secondary MK 
standards in Keenan & McNeil (1989) or have accurate spectral 
types in lKeenan & Pitts! dl98Qh . We find 38 K-type stars and 14 
M-type stars (earlier than M6) classified as giants or supergiants 
in these references. 

3.1. Temperature determination and spectral types 

For late-type stars the dominant features in the /-band spectra are 
the molecular TiO absorption bands, with VO bands also becom- 
ing prominent for spectral types M6 and later (Kirkpat rick et al.l 
1991). These bands are much stronger in luminous stars (gi- 
ants and supergiants) than in stars close to the main sequence. 
The strength of the TiO bands has been used in the past as the 
basis of several temperature calibrations, via the definition of 



diffe rent photometric or spectroscopic indices (see for exam- 
ple lAlvarez & Mennessiedll997l: ICenarro et al.ll2009l) . iRamsevI 
(119811) found that the depth of the TiO 6(0,0) fl-branch band- 
head at 8 660 A correlated strongly with for a sample of K 
and M giants observed at high resolution. Here we build our 
own temperature scale based on the strength of this feature on 
intermediate-resolution spectra. We measure the TiO absorption 
between 8 860 A and 8 872 A against a continuum estimation be- 
tween 8850A and 8 857 A. The magnitude of this jump is then 
calculated as: 



TiO = 1 - 



^TiO 



(1) 



1 The spectra in Cenarro et al. (2001) were observed in different runs 
at several telescopes. Around two thirds of the spectra have resolution 
(FWHM) of 1.2 A, around one third have resolution of 1.5 A, while a 
few have lower resolutions. 



T eff = 3847.56-1 803. 14xTiO 
CT Teff = 1 1 6.36 K 
T eff = 4057.94-9314.81 xTiO 
ct™= 356.06 K 



Cenarro et al. 01 
Rayner et al. 09 




Fig. 3. Calibration of the depth of the TiO bandhead at 8660A 
against T G s for a sample of dwarfs (red circles), giants and su- 
pergiants ( black circles) with M K spectral types present in the 
libra ry of Cenarro et al.l (l200ll) . Adding stars from the library 
of iRayner et al.l (I2009I) . observed at lower resolution (asterisks), 
does not significantly change the trend. 



To calibrate T e s against this measure ment of the TiO band 
strength, we used the spectral library of ICenarro et al.l (2001) 
as a source of primary calibrators. In order to check if there is 
a depen dence on spectral re solution, we also employed the li- 
brary of IRayner et al.l (120091) . This catalogue covers our spectral 
ra nge at R = 2 000, a nd ag ain contains mostly stars classified 
in iKeenan & McNeill dl989b . Both catalogues contain a fairly 
wide sample of spectral types, including giants, supergiants and 
dwarfs. Over the spectral range and at the resolution in question, 
giants and supergiants have very similar spectra. Therefore we 
explored the possibility that a single calibr ation can be used for 
luminous stars of all luminosity classes. As ICenarro et al. I (l200lh 
provide T e ff determinations for all their stars, we used this subset 
to obtain both a relation between TiO band depth and T q q and a 
relation between T e s and G, an index that labels the spectral type 
of the star (G = for a putative O0 star, increasing by one for 
each type, reaching G = 30 for an F0 star and so on). We found 
that for types later than AO, 



log(r eff ) = (4.155 ± 0.008) - (9.79 ± 0.16) x 



1000 



(2) 



Initially, we fit the relationship betwe en the TiO index and 
r e ff for the stars from lCenarro et al.l (12001b . The results are plot- 
ted in Fig. [3] As is seen, the index correlates well with tempera- 
ture for spectral types later than M0; the resultant relation being 



T eff = (3848 ± 38) - (1803 ± 154) x TiO . 



(3) 



Giants and supergiants canno t be separated . 

For stars in the catalogue of IRayner et al.l (120091) . we only 
have spectral types. Consequently, we first estimated their T e s 
via Eq. [2] With this estimate, we then chec ked that they 
follow the same (770,r e ff) relation as stars in ICenarro et al.l 
( 200 l b (Fig. 13 ). Typical errors are ±116 K for the stars from 
ICenarro et al.l (120011) and ±143 K for stars in IRayner et al.l 
( 2009). Both dispersions are fully compatible, as the latter value 
includes the errors from both Eqs. [2] and [3] Giants and super- 
giants follow the same relation, while the difference in spectral 
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resolution between the two samples does not seem to affect the 
index significantlyB 

Main- sequence stars, on the other hand, have rather weaker 
TiO features, and follow a very different relation, 

T eff = (4058 ± 140) - (9315 ± 2190) x TiO . (4) 

Because of the weakness of the TiO band, the relation shows 
an increased dispersi on, with typical errors ±356 K for stars from 
ICenarroetaD (12001b . 



3.2. Luminosity calibration 

We then proceeded to calibrate the main luminosity indicators, 
namely the strength of the Can triplet and the Tii dominated 
blend at 8468A. For this, we measured the E Ws of these fea- 
tures on standard stars in ICenarro et al.l (2001 h using the same 
measuring technique as for our targets. As the number of su- 
pergiants is low in this sa mple, we comp l ement the dataset 
with a few more stars from iGinestet et al.l (1 19941) , at slightly 
higher resolution. The results for the two strongest Ca n lines 
are displayed in Fig. HI We find "normal" luminosity class 
III stars to have EW(Can 8542+Can 8662) in the range 6- 
8 A, with a clear tendency to lower values at early K types. 
Stars whose classification includes metallicity anomalies devi- 
ate slightly from the general trend. Ba stars generally have high 
values, but all giants have EW(Can 8542+Can 8662)<8.5A. 
The few luminosity class II bright giants in the sample have 
EW(Can 8542+Can 8662)^8.5A, while the seven supergiants 
earlier than M3 have EW(Can 8542+ Can 8662)>9A Thes e 
results agree qualitatively with those of IGinestet et al.l (Il994j) . 
though our definition of the continuum leads to smaller EWs0. 
The EW increases on average with spectral type, until it reaches 
a maximum at K4-M1 , and then decreases, due to the depression 
of the continuum lev el by the TiO features The criterion defined 
by iDfaz et al] {l989) to separate supergiants from bright giants 
at solar composition, namely EW(Can 8542+Can 8662)>9A, 
also appears valid for our resolution and continuum definition, 
at least for spectral types earlier than M3. 

Figure [5] shows the measurements of the 8468A blend. We 
find that "normal" class III stars have EW(8468) in the 0.6- 
0.8 A range, again with the lower values corresponding to the 
earlier types. Metallicity anomalies cause important dispersion 
here. In particular, stars classified as CN-strong have values in 
the 1.0- 1.2 A ran ge. These objects are b elieved to be to some 
degree metal-rich (IKeenan'&' Heck 1994). All bright giants have 
EW(8468)=0. 9-1.1 A, with all supergiants displa ying values > 
1.2A. Our values agree very well with those of Ginestet et al. 
(11994 their Fig. 13), as the local continuum in this range coin- 
cides with their selected high points. We can thus confirm that 
the EW of this blend is a very good luminosity indicator, allow- 
ing perhaps a better separation of supergiants and giants of nor- 
mal composition than the Can lines. Interestingly, all the CN- 
strong giants have EW(Can 8542+Can 8662)^ 8 A, high for 



12.5 



2 In order to explore the effect of spe ctral resolution, we a rtificially 
degraded the resolution of the spectra in Cenarro et al. (2001) to lower 
resolutions. We find that the TiO index starts to loose sensitivity to the 
earliest M spectral types at resolving power R ^ 1500 and deviates 
strongly from the relation shown in Eq. [3]for R ^ 1000. 

3 IGinestet et al] d 19941) use a few high points in t he spectra to de- 
fine the cont inuum, while most other references (e.g. JPiaz et alj |l989: 
Mallik||l993) use small featureless sections at the edges of the spectral 
range. As a consequence, Gineste tet al] d 19941) measure higher EWs 
than other authors do. 
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Fig. 4. Measurements of the EW of Can for stars in the li- 
brary of Cenar ro et al. Circles are luminosity class III 
giants (open circles are stars whose spectral classification in- 
cludes some chemical anomaly, while filled circles have no 
such qualifiers). Diamonds are luminosity class II bright giants. 
Triangles are supergiants (triangles pointing upwards have lu- 
minosity class lb, while those pointing downwards have class 
lab). With one single exception, giants and supergiants are well 
separated up to spectral type M2. The strongest Ca n triplet cor- 
responds to the Kl la-lab supergiant HD 63302. 
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Fig. 5. Measure ments of the EW o f Ti i/Fe i blend for stars in 
the library of ICenarro et al.l (12001 1) . Symbols are as in Fig. 5] 
Except for a few deviant points at K3, supergiants are very well 
separated from normal-composition giants. The strongest feature 
is measured again in HD 66302. 



class III, but always lower than supergiantfl Therefore a com- 
bination of EW(Can) and EW(8468) allows the identification of 
metal-rich giants and their separation from supergiants. 

For spectral types later than ~M4, the measurement of this 
feature becomes very difficult, as it is located inside the TiO 
XX 8432-8442-8452 feature. Indeed, Fig. [5] shows high disper- 
sion with luminosity class for stars later than M2, the spectral 
type at which this TiO feature starts to become strong. As a con- 
sequence, the luminosity classification for stars later than M3 



4 This is not surprising, as the high EWs of the Ca n lines in super- 
giants are due to the development of broad wings at low gravity 
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can only be achieved by considering the relative strength of the 
whole metallic spectrum. 

3.3. Radial velocities 

Cross-correlation, either with a model spectrum or a star with 
known radial velocity v ra d, is the most effective way to derive 
radial velocities via the Doppler shift introduced in the spectrum 
by motion in the radial direction. If the template spectrum is 
chosen to match closely the stellar spectrum, accuracies in v ra d 
as high as ^ 10% of the spectral resolution can be obtained for 
spectra including a large number of lines. 

We tried two different approaches to derive radial velocities. 
Initially, we employed direct cross-correlation between our spec- 
tra utilising the hcross task within figaro. This allows direct 
comparison between our targets. Radial velocities may be mea- 
sured relative to members of Ste 2, whose velocities in the Local 
Standard of Rest (LSR) have been accurately measured by lD07l 

Unfortunately, in the spectral range that we have observed, 
the spectrum of late giants and supergiants changes drastically 
with r e ff due to the appearance of strong TiO bands. As a con- 
sequence, we need to use as templates a set of spectra that can 
match our heterogeneous dataset. The stars in Ste 2 do not cover 
the whole range of spectral types of our targets. 

In view of this, we decided to use again the spectral library of 
Cenar ro et al. (2001), taking advantage of the fact that the spec- 
tra cover the whole parameter space in and luminosity at 
a resolution comparable to ours. In addition, the spectra in the 
library have been shifted to a common velocity. 

The algorithm that obtains the Doppler shift is a two-step 
process. It starts with a first estimation of T e ff via the TiO cal- 
ibration presented in Section 13.11 which is then used to select 
a template for comparison from the library. Successive velocity 
shifts, in steps of 1 kms" 1 , are applied to the template, and the 
correlation maximum gives us a first estimation of v ra d. Using 
this value, we shift the observed spectrum to the rest velocity, 
and then a second is derived. This second estimate is in most 
cases within ±1 spectral types from the first one. Although a new 
template could be selected, such a small change in spectral type 
has a negligible effect on the derivation of v ra d- Observed v ra d 
were transformed into the LSR reference system using IRAF's0 
rvcorrect package. 

The measured values of v ra d depend on the spectral range 
chosen. Because of the presence of strong telluric bands, wide 
regions of the spectrum cannot be used in the cross-correlation. 
As the humidity was high during the observations, we restricted 
the spectral range used for cross-correlatio n to 8370-8900 A. I n 
Fig. [21 we display the sky transmission from Hin kle et al.l(l2003l) . 
The strongest atmospheric features in this range have 5% of the 
intensity of the O2 band centred at 8231 A. By scaling against 
the 8231 A band, we estimate that the strongest atmospheric fea- 
tures reach a depth of ^ 0.03 continuum units, while the typ- 
ical metallic lines in the spectrum are close to 0.1 continuum 
units. Therefore no atmospheric features are noticeable in our 
spectra and they should not contribute significantly to the cross- 
correlation function. 

Likewise, interstellar features in this range are not numerous. 
The only strong interstellar feature is the 8 620 A diffuse inter- 
stellar band (DIB), discussed in the next section. There is also a 

5 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation 



weaker DIB at 8763 A (e.g IGredel et al Il200ll) and a broad shal- 
low DIB around 8648 A dMunari et al.l 120081: iNegueruela et al.l 
l2010bl) . Finally, interstellar lines corresponding to the (2,0) band 
of the C 2 Phillips system a re seen in stars with high redden- 
ing (e.g.. IGredel et al.ll200ll) . The strongest features, such as the 
Q(2) and Q(4) transitions at /U8761, 8764 A, reach intensities of 
a few hundredt hs of the continuum in the spectrum of Cyg OB 2 
#12 (A y « 10; IGredel et al1l200lh . but are not seen in any of 
our spectra. This lack of detectability suggests that the cross- 
correlation is not affected by the presence of interstellar features. 

We derived an initial estimate of the accuracy of the mea- 
surements by performing cross-correlations against different 
templates and among our targets. From this we estimated the 
internal error of the procedure to be around ±5 km s" 1 . 



4. Results 

A quick look at the spectra obtained (Figs. [6] to [9]) clearly 
shows that all the stars selected by our photometric criteria are 
luminous late-type objects. Accurate classification can be ac- 
complished using the criteria explored in the previous section 
and direct comparis on to the standard stars in the library of 
ICenarro et al.l (l200ll) . To determine luminosity, we use the EW 
of the Ca 11 triplet lines and the 846 8 A blend, as well as the direct 
comparison of the full meta llic spectrum against s tandard stars in 
ICenarro et aD (1200 ll) and ICarquillat et al.l (119971) . For stars later 
than M3, the determination of luminosity is difficult, and relies 
solely on the apparent strength of the metallic spectrum. 

In addition to spectral classification, two other spectral fea- 
tures are important for the analysis. Firstly, the radial velocity 
of the o bjects can provide in formation about their distance. As 
shown by lDavies et al.l (120081) . the Galactic rotation curve allows 
the separation of stars located in different Galactic arms. In this 
direction, the line of sight cuts the Sagittarius Arm at ~ 1.5 kpc 
from the Sun (where vlsr values are ^ 25 k m s" 1 ) and cro sses the 
Scutum-Crux Arm at a distance -3.5 kpc (lTurner| [l980). where 
vlsr values should be ~ 50 kms -1 . Then it reaches the struc- 
ture close to the edge of the Long Bar, perhaps touching it at 
~5 kpc and crossing it for 2 or 3 kpc. It does not reach any other 
arm until it cuts again the Sagittarius Arm at ^10 kpc. Stars i n 
RSGC1 have an average v LS r = 120 kms -1 (Da vies et al1l2008l) . 
while those in Ste 2 have vlsr =110 km s" 1 (lD07l) . Stars in the 
Scutum Complex thus have much higher radial velocities than 
foreground or background objects. 

Secondly, the DIB centred at 8620A has been found to be 
an excellent tracer of reddening. For moderate reddenings, a 
very tight correla tion E(B - V) = 2 7 2(+0.03) x EW(DIB) (A) 
has been found (Muna ri et al. 2008). The DIB does not show 
any correlation with reddening for stars in the heavily red- 
dened (A v = 10.8) open cluster Westerlund 1 (iNegueruela et al.l 
l2010bh . OB supergiants in this cluster show EW(DIB)= 1.1 ± 
0.1 A, suggesting that the band saturates at this value. Such a 
strong feature is easily recognisable in the spectra of heavily red- 
dened stars and serves as an additional indicator of high distance. 
In RSGs, it appears blended with the Fe 1 8621 A line (see Fig.0. 
Because of this, its EW cannot be measured accurately at our 
resolution. However, the presence of this feature can be deduced 
by comparison to Fei 8611 A. This line is always stronger than 
the Fei 8621 A feature at all spectral types. Observing a ratio 
Fei 8621A/Fei 8611A>1 implies the presence of a moderately 
strong DIB blended with the Fe 1 8621 A feature. 
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Fig. 6. Spectra of seven red supergiants in Ste 2. Note the promi- 
nent DIB at 8620A compared to the neighbouring Fei 8611 A. 
The Fe i 8621 A line can be seen as a bump on the red side of this 
feature, because it is displaced to /18624 due to the high v ra d of 
these objects, except in D3, which is too late to show this feature. 



4.1. Stephenson 2 

W e ob tained spectra of seven stars listed as members of Ste 2 
by lD07l . Their spectra are displayed in Fig. [6] All the objects are 
M-type luminous stars. In Table [2 we list their spectral ty pes as 
estimated from the strength of the CO 2.3/mi bandhead in ID07L 
as estimated using the strength of the TiO 8 660 A band and as 
derived from direct comparison to standard stars. 

The agreement in gen eral i s very good, except for D18, 
which was classified M4 b y lD07l and appears rather earlier here. 
The observations of lD07l were taken in 2006, three years be- 
fore our run. Some RSGs are known to show spectral variability, 
though this is a r are occurrence, perhaps confined to the most 
lu minou s objects (iLevesque et al.|[2005l; IClark et al.ll2010h . 

lD07l find an average radial velocity for the cluster vlsr = 
109.3 ± 0.7 kms" 1 , with the uncertainty representing Poisson 
statistics for 26 members with measured values. The seven stars 
that we have observed have average vlsr = 109+ 2km s" 1 (where 
the error is the standard deviation of the seven measurements), 
i.e., are representative of the cluster velocity. Their individual 
values range from 105 km s" 1 for D29 to 1 12 km s" 1 for D15. 

When we compare our values, we find that the differences in 
v ra d (D07 minus our values) range from + 13 km s" 1 to -2 km s" 1 , 
with an average difference of +8 km s" 1 and a standard deviation 
of 5 kms" 1 . For D3, we have three measurements, which give 



values of 100, 111 & 93 km s \ corresponding to an average of 
101 ±9 kms" 1 . 

Based on these values, we estimate that the total error in one 
single measurement of v ra d using our technique is ±10 km s" 1 or 
slightly higher^. However, the technique appears robust when 
applied to a sample of objects in a cluster. We assume that the 
difference of +8 kms -1 represents our systematic shift with re- 
spect to the LSR system adopted by ID07L based on the radial 
velocity of Arcturus. The average radial velocity that we obtain 
for Ste 2 is vls r = 102 ± 4 kms -1 . The difference with the av- 
erage value of |P07l is +7 kms -1 , consistent with a systematic 
shift. 

The reason for this shift is unclear. It may be due to a 
different definition of the LSR. rvcorrect transforms using the 
Standard Solar Motion (+20 kms -1 towards I = 56°, b = 23°). 
We have calculated mean velocity corrections for the stars ob- 
served in Ste 2 AV© = -4.5 kms" 1 and AV LS r = +16.0 kms" 1 . 
In any case, the shift is of the same order of magnitude as 
the dispersion in our measurements for cluster members, and 
so could be statistical, as the number of members observed is 
low. In addition to the accuracy of individual measurements, 
a real dispersion in v ra d amongst cluster members could be 
present. Observations of five RSGs in the association surround- 
ing x Persei show random fl uctuations in rad ial velocity with 
amplitudes > 5 kms" 1 dMermilliod et al.ll2008l) . 

4.2. RSGC3 

We observ e d five stars in RSGC3, those numbered Sl-5 in 
IClark et aTl (l2009al) . Their spectra are displayed in Fig. [7] The 
five objects are M-type supergiants. Their measured v ra d are con- 
sistent within the accuracy discussed above (S5 shows the value 
deviating most clearly, but it also has the poorest spectrum with 
SNR-35 per resolution element). The average velocity for these 
five members is vlsr = +95+7km s" 1 , where the error repr esents 
the standard deviation. Interestingly J Alexander et al1(l2009l) find 
that the CO distribution in this direction shows a maximum at 
vlsr = +95 kms" 1 . Material with this vlsr shows a spatial dis- 
tribution strongly correlated with the GLIMPSE 8 jim emission, 
and presents a hole around the cluster, suggesting a cavity blown 
by stellar winds and supernova explosions. The vlsr for RSGC3 
is smaller than that found for Ste 2 (vlsr = 102 + 4 kms" 1 ), 
though consistent within the errors. 

With these new spectral types, we calculate in dividual colour 
excess es E(J - K$), using the calibrations of lLevesque~ et al. 
( 2005|) to calculate the intrinsic (/ - K)o f° r each spectral typeQ. 
Th e values der i ved, sh own in Table [3] are very similar to those 
of [Clark et al.1 (l2009al) . Given the relatively good agreement 
between our spectral types and those derived using the CO- 
bandhead calibration, changes in the individual absolute magni- 

6 This value is slightly higher than expected. However, we note that 
it is not due to use of external templates, as direct crosses between the 
three spectra of D3 using hcross give similar dispersion. We can only 
speculate about the reason why the error is larger than the expected 10- 
15% of the spectral resolution. It may be related to the fact that we tried 
to align several stars in each slit, meaning that objects are occasionally 
not well centred, perhaps resulting in different light paths for starlight 
and the arc lamp light. 

7 Note that the (/ - K) c olours obtained by applying the calibration 
of r eff versus spectral type in lLevesque et al.l (120051) to their (J-K^/T^ 
relationship are almost identical to the empirical calib ration of intrinsic 
(J-K) for red giants of Str aizvs & La zauskaite (2009). This agreement 
not only gives support to both calibrations, but indicates that the (J-K) 
colours of red giants and supergiants are indistinguishable. 
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Table 2. New spectral types for members of Stephenson 2. 



ID 


Co-ordinates 


i 


J 


H 




Spectral 


Spectral 


Spectral 


Vlsr 




RA 


Decln. 










Type (CO) 


Type (TiO) 


Type 


/t — 1 \ 

(km s l ) 


D3 


18 39 24.6 


-06 02 13.8 


12.85 


7.27 


5.46 


4.50 


M4I 


M6 


M5 lab 


+ 101 


Dll 


18 39 18.3 


-06 02 14.3 


14.35 


8.35 


6.21 


5.26 


M4I 


M3 


M3.5Ib 


+ 106 


D15 


18 39 22.4 


-06 01 50.1 


13.07 


8.13 


6.35 


5.51 


M2I 


M2 


Ml. 5 lab 


+99 


D16 


18 39 24.0 


-06 03 07.3 


13.20 


8.24 


6.44 


5.60 


M3I 


Ml 


Ml la 


+97 


D18 


18 39 22.5 


-06 00 08.4 


12.88 


8.18 


6.45 


5.63 


M4I 


Ml 


Ml. 5 lab 


+ 101 


D20 


18 39 24.1 


-06 00 22.8 


13.18 


8.43 


6.66 


5.81 


M2I 


Ml 


Ml. 5 lab 


+ 103 


D29 


18 39 22.2 


-06 02 14.7 


13.17 a 


8.61 


6.88 


6.15 


M0 1 


Ml 


M0 lab 


+ 107 



{a) The DENIS magnitude seems to correspond to a blend of two stars, D29 and 2MASS J18392193-0602165. 



Table 3. Properties of confirmed members of RSGC3. 



ID 


Co-ordinates 


r 


J a 


H a 


K s a 


Spectral 


Spectral 


Spectral 


E(J-K S ) 


Vlsr 




RA 


Decln. 










Type (CO) 


Type (TiO) 


Type 




(kms" 1 ) 


SI 


18 45 23.60 


-03 24 13.9 


13.5 


8.55 


6.54 


5.58 


RSG 


M3 


M3.5Ib 


1.80 


+97 


S2 


18 45 26.54 


-03 23 35.3 


b 


8.53 


6.62 


5.75 


M3I 


M2 


M2Ia 


1.68 


+ 104 


S3 


18 45 24.34 


-03 22 42.1 


13.7 C 


8.54 


6.43 


5.35 


M4I 


M4 


M4I 


2.02 


+99 


S4 


18 45 25.31 


-03 23 01.1 


13.7 C 


8.42 


6.39 
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(a) JHK S magnitudes from 2MASS. / magnitudes from USNO-B1.0. 

No USNO-B1.0 source, but the star is similarly bright to the other four members observed in the DSS2 / image. 
(c) The magnitude corresponds to a blend of two stars in the DSS2 / image. 



4.3. Field 1 

In Field 1, we selected 12 candidates (listed in Table[T]) and man- 
aged to observe seven of them. All the stars observed turn out 
to be M-type luminous stars. Two objects, F1S01 and F1S10, 
have relatively low E(J - Ks) and show no evidence for heavy 
extinction in their spectra, as the 8 620 A DIB is not prominent. 
Both are quite late and their exact luminosity classes are dif- 
ficult to assign, but both look rather luminous, close to being 
low-luminosity supergiants. 

The other five objects have similar values of E(J - Ks), all 
within the same range as stars in RSGC3. As seen in Fig. [8] all 
five show very prominent 8620A DIB, with EW^1.3A. The fea- 
ture is so broad that it is almost blended with Fei 8611 A. All 
five objects have high observed v ra d, but our measurements show 
higher dispersion than those for stars in Ste 2 and RSGC3. The 
average of their velocities is vlsr = +89 + 16 kms" 1 , where 
the error represents again the standard deviation. All the individ- 
ual measurements are within one standard deviation, but deviate 
slightly more than the estimated error of + 10 - 12 km s" 1 (all the 
spectra are moderately poor with typical SNR-50). 

Among the objects that we did not observe, F1S03 with 
(J - K s ) = 1.89 and FISH with (/ - ^ s ) = 1.79 must have rel- 
atively low colour excesses if they are K-M stars. On the other 
hand, F1S05 and F1S09 have colours similar to the group of five 
supergiants. Below we argue that all these high-reddening ob- 
jects represent a new cluster of red supergiants. 



Wavelength (A) 4.4. Field 2 

Fig. 7. Spectra of five red supergiants in RSGC3 . Note the promi- xhere is an obvious concentration of bright stars in the K s image 
nent DIB at 8620A compared to the neighbouring Fe i 861 1 A. f the area. Only 1 1 of them pass our photometric cut, defining 

an apparent compact group around (RA: 18:44:45, 6\ -03: 14:30) 
and a smaller group ~ 6' to the North. We find seven candidates 
tudes are small. Errors in M K are dominated by the uncertainty in the main group and four in the northern group. Their coordi- 
in the distance to the cluster. nates and magnitudes are listed in Table [T] 
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Fig. 8. Spectra of red supergiants in Field 1. The five top ob- 
jects have similar colours and are affected by heavy reddening. 
The bottom spectrum corresponds to S10, which is less lumi- 
nous than the other objects. It has much lower reddening than 
the other stars (as shown by the weak 8620A DIB) and very dif- 
ferent radial velocity. 



We observed five candidates in the main group. Of these, 
three objects F2S01, F2S03 and F2S04 are M-type supergiants 
with similar colours. Their colour excesses are almost identical, 
and similar to the lowest values displayed by stars in RSGC3. 
Their observed v ra d are also similarly high. However, there is no 
strong evidence for a large cluster here. F2S05 and F2S07 are 
both late-M giants. Their colour excesses are much lower, even 
if their v ra d are high. Of the stars that we did not observe, F2S06 
has (J - Ks) = 1.56, and cannot be highly reddened if it is a red 
star. F2S03, on the other hand, has (J - Ks) = 3.46, and it could 
be a redder or more highly reddened supergiant. 

In the northern group, F2S09 has a late spectral type. The 
spectrum is quite noisy, and the luminosity class is difficult to 
determine. It looks like a giant (see Fig. [9]), but displays high red- 
dening and high v ra d, similarly to the supergiants. F2S11, which 
we did not observe, has colours and magnitudes similar to the 
three supergiants in the main group. 

4.5. Field 3 

We observed six of the ten candidates selected in this area 
(Table []]). All of them are mid or late M giants, with moder- 
ate reddening. Though their colour excesses E(J - Ks) are quite 
uniform, they display a very wide range of observed v ra d- 




8400 8500 8600 8700 8800 8900 

Wavelength (A) 



Fig. 9. Spectra of stars in Field 2. The three top objects have sim- 
ilar colours and are affected by heavy reddening. They are likely 
members of the RSGC3 association. F2S05 is a luminous fore- 
ground giant. F2S09 morphologically looks like a giant (though 
the spectrum has low SNR), but its reddening and v ra d are similar 
to those of supergiants in the RSGC3 association. 

The four stars that we did not observe have redder colours 
than these giants, with F3S04 and F3S07 showing 2MASS 
colours similar to members of RSGC3. 



5. Discussion 

We have observed a sample of RSG members of the massive 
clusters Ste 2 and RSGC3 in the region of the Can triplet. In 
addition, we have observed 19 objects within 20' of RSGC3, 
selected because their 2MASS colours suggested that they might 
be reddened luminous late-type stars. 

5.1. Validity of criteria and calibrations 

We have selected targets based on their infrared brightness and 
colours. The criterion 0.1 < Q < 0.4 was expected to reject 
early-type stars and red dwarfs. Indeed, all the 19 objects se- 
lected according to these criteria for which we have spectra turn 
out to be luminous M-type stars with different amounts of red- 
dening. Eight objects are highly reddened supergiants. The rest 
are mid and late M giants or low-luminosity supergiants with 
different amounts of reddening, ranging from E(J - Ks) = 0.55 
for F1S01 to 1.36 for F2S09. Therefore we may conclude that 
the criteria defined in Section |2J"1 to identify reddened luminous 
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late-type stars are very successful, as all the objects selected are 
M supergiants or late M giants. The selection in Q effectively re- 
moves most (perhaps all) low-mass red stars, which are expected 
to have Q > 0.4. The spectral characteristics of the objects se- 
lected suggest that their lower Q values (lying in the region cor- 
responding to yellow stars) are likely due to the existence of ex- 
tended atmospheres, which modify their spectral energy distri- 
bution. As the criteria used are mainly related to intrinsic stellar 
properties rather than to the characteristics of the area surveyed, 
these criteria may be effective at selecting obscured red luminous 
stars in other regions. 

Unfortunately, the photometric properties of luminous red 
giants and RSGs seem to be identical, and a separation of both 
classes is not possible based solely on 2MASS photometry. With 
the small sample observed here, it is not possible to assess if a 
combination of 2MASS and Spitzer photometry can effectively 
separate the two classes. Though our sample shows some evi- 
dence that this may be the case, objects enshrouded in dusty en- 
velopes will necessarily display very red colours caused by dust 
emission, independently of their luminosity. A larger sample of 
objects with accurate spectral types will thus be needed to assess 
this issue. 

A second important result concerns the use of the CO band- 
head features in the K band to estimate the spectral type. 
Different authors have calibrated the EW of the bandhead against 
spectral type, finding a degree of confusion between giants and 
supergiants, and an e stimated uncertainty of +2 spectral sub- 
types for a given star (iDavies et al.ll2QQ7l I Alexander et al.l l2009: 
Negueru ela et al.ll2010al) . Our targets include 1 1 red supergiants 
in Ste 2 and RSGC3 for which previous authors had estimated 
the spectral type using the calibration of the CO bandhead. If we 
compare the spectral types estimated with those derived by di- 
rect comparison to the standards, the agreement is excellent. The 
average difference is one spectral type. Since the classification 
of the standard stars was originally carried out in the blue- violet 
region (meaning that the assignment of spectral types based on 
Z band spectra is indirect), the agreement is surprisingly good. 
Even though a given star may have its spectral type misestimated 
by two subtypes, there does not seem to be any systematic ef- 
fects. 

Finally, for the 3 1 stars in our sample, we can also compare 
the spectral types de rived from the calibration of the TiO band 
described in Section 13.11 with those derived by direct compari- 
son to the standards. In this case, the average difference is only 
0.6 spectral types, confirming that this calibration provides valid 
spectral typing for M-type giants and supergiants. 

5.2. A new cluster: Alicante 7 = RSGC5 

Since five RSGs detected in Field 1 seem to have very similar 
colours, we investigated the possibility that they comprise a new 
cluster of RSGs. To accomplish this, we unde rtook a similar pho- 
tometric analysis similar to that presented in (iClark et al.l l 2009a) 
for RSGC3. We took a 7' circle around a nominal centre, chosen 
as the position of F1S06 and selected 2MASS_photometry for 
stars with quality flags "AAA" and AK S < 0.050. 

We found a number of bright objects, clearly grouped in the 
(/ - Ks)/K$ and (H - K$)/Ks diagrams, which stand up from 
the rest of the field (see Fig. [TO]). This group includes all the five 
M-type supergiants identified, plus F1S09 and three other stars 



4 - 



8 Note that S8 does not fulfil these conditions, as it is saturated in 
K s . Nevertheless, we include it in the sample, as it is certainly a red 
supergiant. 
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Fig. 10. 2MASS colour magnitude (/ - ^ s )/^s plot for stars 
within 7' of the centre of the new cluster Alicante 7. The dia- 
monds identify stars with colours and Q in the same range as 
RSGs in RSGC3. Five of them, marked by crosses, are spectro- 
scopically confirmed as M-type supergiants. Triangles are other 
stars in the area with very similar colours, but fainter magni- 
tudes. Note that the brightest star, A6 in Table [4] is saturated in 
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Fig. 11. Finding chart for Alicante 7, with the stars listed in 
Table |4] indicated. The finder comprises a K-band image from 
2MASS with a ~ 11' x 11' field of view centred on the po- 
sition of A5 (=F1S06 in Table [T]). Red circles indicate likely 
RSG members (including the five objects spectroscopically con- 
firmed), while blue squares mark stars with similar colours, but 
fainter magnitudes. 



with similar colours and the expected Q, but faint / magnitudes. 
In addition, we found four other stars with similar colours but 
fainter magnitudes. This second group is not so clearly sepa- 
rated from the field population as the nine main candidates. We 
identify these nine main candidates as a new cluster of red su- 
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pergiants, which we provisionally name Alicante 7. Figure Qj] 
shows a finder for all the candidates. 

Star A6 (=F1S08) is the brightest cluster member in the 
band. Its magnitude is in the saturation range for 2MASS. 
This object is a strong mid-IR source, and almost certainly the 
counterpart to IRAS 18418-0331. This suggests that it may be 
surrounded by a dusty envelope. If this is the case, its luminosity 
class, based on the strength of the metallic spectrum, may be 
incorrect, as the spectrum may be veiled by dust, leading to the 
dilution of the metallic features. 

We can obtain an estimate of the distance to Alicante 7 by 
studying the distribution of inter stellar extinction in it s direc - 
tion. Following the technique of llCabrera-Lavers et aD (120051) . 
we used the population of red clump giants in infrared 
colour/magnitude diagrams. As the 2MASS data are not deep 
enough for this purpose, we used UKIDSS data. We plotted the 
(/ - K)/K diagram within 20' of the cluster centre, and selected 
the giant population, which is seen in such plots as a well defined 
strip. We chose a 20' radius as a compromise that allowed us 
to construct a representative (d, A K ) curve for the cluster sight- 
line, while providing a large enough number of stars to permit 
an accurate determination. Decreasing this value to, for exam- 
ple, l(y produces noisier results, but does not change the overall 
behaviour of the extinction. We perform the same calculation for 
the area centred on RSGC3. Given the distance between the two 
clusters, the two areas partly overlap. The results are shown in 
Fig. El 

As expected, both extinction runs are similar. The extinction 
starts at low values and suffers an important increase at a dis- 
tance of 3.5-4 kpc, which is believed to correspond to the posi- 
tion of the first crossing of t he Scutum- C rux Arm along this line 
of sight. Around £ - 28°. lTurner](ll980l) found several luminous 
OB supergiants with distances in excess of 3 kpc and reddenings 
E(B - V) « 1.3, corresponding to A K - 0.4, in very good agree- 
ment with our determination. After this, the reddening increases 
steadily, reaching a value A K « 1.5 at 6 kpc. The extinction in- 
creases abruptly between 6 and 7 kpc, reaching values that ren- 
der the technique unusable for higher distances. From this, it is 
clear that both Alicante 7 and RSGC3, with A K « 1.5 are lo- 
cated at a distance ~ 6 kpc, just in front of the extinction wall. 
This distance determination i s in very good agreem ent with the 
dynamical distance found by I Alexander et al. I (12009) for the CO 
material apparently surrounding RSGC3, namely d = 6.1 kpc. 

5.3. An extended association surrounding RSGC3 

Comparison of Tables [3] and @] reveals that the RSGs in RSGC3 
and Alicante 7 have very similar colours and magnitudes, sug- 
gesting that they share a common distance and age. In Fig. [13] 
we plot the 2MASS (/ - ^s)/^s diagram for all the stars with 
high-quality photometry within 7' of the centres of both clus- 
ters, together with all our photometric candidates from TableQ] 
Several populations can be seen in both diagrams, revealing the 
similarity of both sightlines. At 6 kpc, the separation between 
the two clusters is 24 pc. Therefore the two groups of RSGs must 
be two separate clusters. Both clusters, however, appear behind 
approximately the same amount of obscuration, and their RSGs 
display very similar 2MASS colours and magnitudes. Given that 
the two clusters are situated at the same distance, the RSGs must 
have the same intrinsic luminosity, indicating that the two clus- 
ters have the same age. 

The age of R SGC3 has been estimated at 16-20 Myr 
dClark et al.ll2009al) . Simulations of stellar populations using a 
Salpeter IMF suggest that a population of 10 000 Mq at this age 
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Fig. 13. Composite 2MASS colour magnitude (/ - ^ s )/^s plot 
for the whole area surrounding RSGC3. Black dots are stars 
within 7' of the centre of RSGC3, while blue microdots are stars 
within 7' of the centre of Alicante 7. Red squares represent the 
16 most likely RSG members of RSGC3, most of which are con- 
firmed spectroscopically. Yellow diamonds as likely RSG mem- 
bers of Alicante 7, as in Fig. [TO] Blue triangles are all the stars 
selected in Field 2 and Field 3 with i > 13.2. The three ob- 
jects marked with crosses are the three M-type supergiants with 
high reddening spectroscopically identified in Field 2. Finally, 
the circles represent all the objects selected in the three fields 
with i < 13.2, none of which is a reddened supergiant. 

should contain eight RSGs dClark et al.ll2009ah . Other simula- 
tions, using a Kro upa IMF, reduce th e number to 2-5 RSGs for 
each 10000 M Q dClark et al.ll2009bl) . Based on this, Alicante 7 
must contain a minimum of 10 000 M© to have a population of 
nine RSGs and could be up to ~ 3 times more massive. 

The small group of RSGs found in Field 2 may be a smaller 
cluster, with four or five supergiants, but this is difficult to 
demonstrate; we refer to this putative cluster as Alicante 9. As 
long as the unevolved populations of the clusters in the Scutum 
Complex cannot be detected, it will be impossible to distinguish 
similarly small clusters from random concentrations of bright 
red stars. As a consequence, estimating the total mass of the as- 
sociation surrounding RSGC3 is very difficult. For example, a 
cluster with only four RSGs, which will be very difficult to iden- 
tify as such, may easily contain ^ 5 000 Mq. 

In order to make a rough estimate of the population, we con- 
sider the population of bright red stars that are not members of 
the association. Of the 19 objects in the area surrounding RSGC3 
for which we have spectra, 1 1 are unlikely to be members, as 
they present much lower extinction than cluster members. Most 
of them have E(J - Ks) ~ 0.6-0.8 and relatively low vlsr- 
They are likely located at distances d ~ 3.5 kpc, where the line 
of sight cuts the Scutum-Crux arm for the first time0. All these 
objects are relatively bright in i, with typical values in the 1 1-12 
range, in good agreement with the moderately low extinction^. 

9 At this distance, t he Galactic rotatio n curve predicts velocities 
vlsr ~ 40 - 60 km s _1 dPavies et alj |2008\ meaning that there is good 
agreement between reddening and v LS r for a number of objects such as 
F1S01, F1S10, or F3S08. Other objects, such as F3S01 or F3S05, have 
very discrepant v LS r 

10 One exception is F2S09, which has high E(J - K s ) and v LS r = 
94 km s _1 . This object is also the faintest amongst the giants in /.These 
values suggest that it must be just in front of RSGC3. 
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Fi g. 12. Left panel: Run of th e extinction in the direction to Alicante 7. The data have been obtained by applying the technique 
of ICabrera-Lavers et al.l (120051) to the red clump giant population within 20' from the nominal cluster centre. Right panel: Same 
for the field surrounding RSGC3. Both runs are pretty similar. The sudden increase in the reddening at ~ 6 kpc provides a direct 
estimate on the distance to both clusters. 



In contrast, members of RSGC3 and other M-type supergiants 
in Alicante 7 and Alicante 9 have rather fainter i magnitudes, in 
most cases ^ 13.5. 

In view of this, it seems possible to separate members of 
the RSGC3 association from interlopers using their i magni- 
tudes. Based on the observed magnitudes of confirmed members 
and confirmed foreground stars, we may set the dividing line 
at i = 13.2. In Fig. [T3} we plot all the photometric candidates 
from the three fields, using different symbols for stars brighter 
or fainter than i > 13.2. We find that stars with i > 13.2 fall 
inside the area of the diagram occupied by cluster members and 
confirmed supergiants. In contrast, stars with i < 13.2 occupy a 
different region of the diagram, with lower (J - K) colours. The 
only exception is F2S09, whose nature is unclear. It has high 
E(J - Ks) and vlsr = 94 kms -1 . Therefore, if it is a giant, it 
must be located just in front of the RSGC3 association. 

Since brightness in i and (J-K$) seem to be well correlated, 
it seems safe to assume that any star in the vicinity of RSGC3 
with 0.1 < Q < 0.4 and 2.5 < (/ - K s ) < 3.7 has a very high 
probability to be a RSG in the RSGC3 association If we turn 
now to the candidate members of RSGC3 listed by IClark et aD 
(2009a), we find that 1 8 of them fulfil all criteria, defining a com- 
pact core ~ 3' across and a cluster radius ~ 6'. 

At larger distances, apart from the nine candidate members 
of Alicante 7, there are ^ 25 objects fulfilling these criteria 
within r = 18' (3 cluster radii) of the cluster core . None of 
them can be identified with a DENIS or USNO-B1.0 sources at 
i < 13.2, but 18 have counterparts fainter than this limit (red 
circles in Fig. []]). Most of these objects are likely RSGs in the 
association. If we take into account that criteria are not designed 
to select all RSGs and that some objects pass all the 2MASS cri- 
teria but have no optical counterpart, the association including 
RSGC3 must contain <: 50 RSGs within 1 8' (~ 31 pc at 6 kpc). 
Using the estimates of IClark et al. (2009b), such a population of 
RSGs would require an underlying population of at least 10 5 Mq. 
Of course, such a huge structure needs to be confirmed via high 
resolution spectroscopy, which can provide kinematic member- 
ships. 

An interesting parallel is provided by the G305 star-forming 
region. Two massive clusters (Danks 1 & 2), separated by only 



3 pc, have very similar, but measurably different ages, and are 
surrounded by a halo of active star formation extending up to 
distances ^30 pc away (IClark & Porterll2004l IClark et aljboilb . 
Such a region may evolve over the next ~15 Myr to resemble 
the association around RSGC3 (though with the mass scaled 
down by perhaps a factor 4-5). On an even smaller scale, the 
Cas OB 8 association, in the Perseus Arm, provides a better 
match for the age of the RSGC3 association. The largest clus- 
ter in the area, NGC 663, is relatively massive, containing (at 
an age 20-25 Myr) five blue s upergiants in the core region and 
at lea s t two RSGs in the halo dMarco et al.l 120071: iPandev et aD 
120051: iMermilliod eTaD l2008h . The smaller cluster NGC 654, 
which contains one blue and one yellow sup ergiant, is located 
~ 30 pc away (for an assumed d = 2.4 kpc: lMarco et al.ll2007b 
IPandev et al.ll2005l) . An even smaller cluster, NGC 659, which 
does not contain supergiants, is -23 pc away from NGC 663, 
while NGC 581, which contains one blue and one red supergiant, 
is ~ 65 pc away. All the clusters in the area have approximately 
the same age, and there seems to be a number of s upergiants at 
a com parable distance spread between the clusters (Hump hreysl 
fl978h . 

6. Conclusions 

We have conducted a pilot study of the neighbourhood of 
RSGC3. We have proposed a number of criteria that allow the 
separation of candidate luminous red stars from 2MASS pho- 
tometry and observed a sample of the candidates in the 8300- 
9000 A range at intermediate resolution. All the objects spectro- 
scopically observed were found to be M-type giants or super- 
giants, confirming the validity of the photometric criteria. Our 
main results are as follows: 

- We confirm the validity of the different calibrations based on 
the strength of the CO bandhead in the K band to estimate 
spectral types of RSGs. The accuracy of ±2 spectral types 
estimated by lD07l seems a good estimate. 

- We confirm membership of five RSG candidates in the open 
cluster RSGC3, and find that their average vlsr is compatible 
within the uncertainty of our method with the average vlsr of 
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Ste 2. Our average velocity for RSGC3, v ra d = 95 + 7 km s" 1 , 
when compared to the rotation curve of the Galaxy in this 
direction, indicates that the cluster is located at d ^ 5 kpc. 

- We find several RSGs in the immediate vicinity of RSGC3. 
Eight of them have colour excesses and vlsr compatible with 
those of cluster members. Photometric criteria suggest an 
even higher number of RSGs in the surrounding area, with 
<: 50 supergiants between the clusters and the association. 

- Five of the RSGs are strongly concentrated and define a can- 
didate cluster, which we call Alicante 7. Analysis of 2MASS 
photometry for this cluster suggests that it contains at least 
nine RSGs. Analysis of the run of extinction in the direction 
to Alicante 7 and RSGC3 shows the distribution of absorbing 
material to be very similar in both directions, and indicates 
that both clusters are located at a distance d « 6 kpc. The 
RSGs in both clusters have identical colours and luminosi- 
ties, suggesting th at they have the same age r=16-20Myr 
(IClark et al.ll2009ah . 

The data presented here strongly point at the existence 
of an extended massive association surrounding RSGC3. If 
higher-precision dynamical data can confirm the association of 
RSGC3 and Alicante 7, together with other objects around the 
clusters, this structure would contain ^60000 Mq, and likely 
^100 000 M Q . 

Our results offer good prospects for a wide-area search us- 
ing high-multiplexing spectrographs. The techniques explored 
in this work offer an effective way of selecting interesting can- 
didates and adequate criteria for the calibration of both spectral 
type and luminosity, which can be easily automated. In a future 
paper, we will present the results of such a survey. 
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